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We studied different levels of resummations of the exactly solvable Bloch-Nordsieck model in order 
to be able to compare the approximations with an exact result. We studied one-loop perturbation 
theory, 2PI resummation and Schwinger-Dyson equations truncated in a way to maintain Ward- 
identities. At all levels we carefully performed renormalization. We found that although the 2PI 
resummation does not exhibit infrared (IR) sensitivity at the mass shell, as the one-loop perturbation 
theory does, but it is still far from the exact solution. The method of truncated Schwinger-Dyson 
equations, however, is exact in this model, so it provides a new way of solving the Bloch-Nordsieck 
model. This method can also be generalized to other, more complicated theories. 

T— I ■ 

o 

(N 
o 

D . I. INTRODUCTION 
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\ In field theories we often encounter infrared (IR) divergences of different kind. Some of them have physical meaning 

■ while others are just artifacts of the perturbation theory. A common feature, however, is that the IR divergences 
' appear at each order of perturbation theory, usually being more and more severe at higher loop orders. In order to 

r~| . consistently define the theory, these IR divergences must be summed up. 

2PI resummations provide a consistent resummation framework known for a long time The basic idea is to 
O ^ replace the free propagator in the perturbation theory with the exact one which is approximated self-consistently with 
D fixed-loop skeleton diagrams. The so-defined perturbation theory is renormalizable and can be applied to 

^ study different physical questions from non-equilibrium [1] , , thermodynamics @ , [I0j|.[l4|.[l5| and different systems 
like 0(N) model [ij] or gauge theories [13|. 

Although the 2PI approximation is constructed by physical arguments, and we expect better results (ie. closer 
\ to the exact one) after 2PI resummation, a priory it is not sure that one really achieves this goal. Probably the 

■ finite lifetime effects are well represented by 2PI resummation both in equilibrium as well in non-equilibrium, 
where the 2PI is close to the Boltzmann-equation approximation [31 . But if the deep IR regime is important where 

^>fi^ multi-particle excitations also play crucial role, the picture is far to be so clean. To make the case even worse, in most 
gauge theory models there is hard to make exact statements about the IR behavior of the model. 

In this paper we aim to study the Bloch-Nordsieck model, which is an exactly solvable 3-l-lD gauge theory \vf\ . 
It is the eikonal approximation model of the QED, and one can argue 19] that in the deep IR regime it describes 
correctly QED. Therefore it is used to estimate IR properties of QED from this model, for example the soft photon 
creation processes [l^l or finite temperature electron damping 

This model is therefore a motivated case study where the accuracy of the 2PI resummation can be assessed. We 
. , therefore perform a series of approximations for the Bloch-Nordsieck model: a one-loop perturbation theory, a 2PI 

■ resummation and finally the solution of the Schwinger-Dyson equations with an Ansatz for the vertex function. In 
this model all levels of the approximations can be treated analytically. We show that the last method is exact in 
the model under consideration - although that is not expected in general. This observation, however, leads us to a 
proposal how the 2PI resummation can be improved in a generic model in order to catch the IR physics correctly. 

The structure of the paper is as follows. We first overview the standard solution of the Bloch-Nordsicck propagator 
at zero temperature in Section |TT1 Then we compute the one loop level fermion propagator in Section IIIIl Next, in 
Section IIVI we write up the 2PI equations for this model, perform renormalization on that, and give the analytical 
solution as far it can be done. Then we study the problem also numerically, determine the fermion propagator and 
compare the result with the exact one. Finally, in Section |V] we study the truncated Schwinger-Dyson equations, 
renormalize it, and show that for the Bloch-Nordsieck model the so-defined approximation scheme is exact. For 
conclusion (Section IVip we summarize our approach again and speculate about the possible generalizations. 
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II. THE MODEL AND THE EXACT SOLUTION 



The Bloch-Nordsieck model is a simplification of the QED Lagrangian where the Dirac matrices 7'^ are replaced by 
a four-vector 

C^^^Ff^.F^^" + ^{iUf,D^-m,)^, iD^=id^-eA^, F^, ^ d^A, - d,A^. (1) 

The singled-out four vector Ufi represents the velocity of the rest frame of the fermion. The fermion wave function 
here has only one component and ^ = ^P*. 

We are interested in the fermion propagator which reads in the path integral representation 

ig{x) = (r*(x)§(0)) = ^ j P§OTPA^e'^[*'*^^l*(x)4'(0). (2) 

At the tree level it reads: 

Qo{p) = , ^ ^ . ■ (3) 

Since it has a single pole, there is no antiparticles in the model, and also the Feynman propagator is the same as the 
retarded propagator. The lack of antiparticles also means that all closed fermion loops are zercO- As a consequence 
the photon self-energy is zero, the free photon propagator is the exact one. In Feynman gauge therefore the exact 
photon propagator is 

G,Ak) = T^- (4) 

Now we shortly review the exact solution for the fermion propagator, cf. [13], [3. We first define the gauge field 
dependent propagator: 

g{x-A) = j 2?^2?*e''5[*'*"^l*(x)§(0). (5) 

This satisfies the Schwinger-Dyson equation 

(m^a'^ - eu^A^ ~ m)g{x- A) = -5{x). (6) 

We will need the solution in an exponential form for A, and this is achieved by introducing the auxiliary quantity 
U{x, v) which satisfies 

i^^^^^^^iiu^.d^' -eu^A^" ^m)U(x,v), U(x,Q)^5{x). (7) 

By integration of the above equation by v and assuming U{x,oo) = (for which we need a convergence factor ie) we 
see that 

g{x) = ~i j dvU{x,v). (8) 

We perform Fourier transformation and separate the free time dependence U(p,v) = e~^'^'^i^P''~"^^'^U{p,v), then we 
obtain 



dU{p,v) f d^k 



leu^J ^A^ip~k)e'^''^P--''-^Uik,iy), f7(p,0) = l. (9) 



The linear — /c^ behavior is the consequence of the linearity in the kernel. If the kernel is non-linear or it is not 
scalar (has a matrix structure) then this form is not true any more. From here an inverse Fourier transformation 
yields 

dfl(r I/) - - }dv'A(uv') 

' =ieu^,A^'{x + uu)U{x,v) => U{x,i^) ^ e° <5(x). (10) 



^ This statement can be best seen in real time representation. There a chain of fermion propagators, because of the retardation, is 
proportional to 0(ti —t2)... ©{tn-i — tn)- In a closed loop tn = ti, therefore the product of theta functions is zero almost everywhere. 
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Once we have an exponential representation for the background-dependent propagator, we can perform the Gaussian 
A-integration. As a result we obtain in the exponent the factor 



2 



- J ^ R:{k)G^Ak)RAk), R^k) = I du'e-"'-^'-'' . (11) 





This integral is UV divergent; in dimensional regularization one finds the result 

a a 



\niuu, (12) 

27r£ TT ^ ^ ' 



where a = e^/ (47r) and /x = •\/47re'*'^/z. Then the fermion propagator reads 



g(p) = -ie^ / rf,^e-'''("''P''-™)+-^"*'^'' = ^-TTX^ (13) 







where Z = T{1 + ^)e2"e/i° . This is UV divergent which means that we need a wave function renormalization. The 
renormalized propagator reads 



{unP" - m)^ 



where C is a finite quantity. 

We can determine the discontinuity of this formula, for simplicity choosing it^ = (1, 0, 0, 0): 

q{p) = BiscGip) = e{po - m) (15) 

oo 

With this spectral function the sum rule / dpog{po) = 1, which is the consequence of the equal time anticommutation 

— CxD 

relations, cannot be fulfilled, since the integral is divergent. This divergence should be compensated with the choice 
C = 0, but then we are faced with a x oo expression. Therefore one should always use a regularized version of the 
spectral function (or propagator), maintaining the sum rule, and only at the end of the calculation is one allowed to 
release the regularization. 

The Lagrangian is Lorentz-invariant in the sense that we must also transform u. So we can choose a Lorentz- 
transformation where Au = (mq, 0,0,0). If is a 4-vclocity then uq = 1; if it is of the form u = (l,v), then it is 
uo = \/l — v^. After rescaling the field — >■ / y/uo and the mass as m — >■ uom, the Lagrangian reads 

£= -jF^,F'''^ + *(i£>o-m)*. (16) 
This Lagrangian will be used mostly in this work later. If necessary, the complete u dependence is easily recoverable. 

III. ONE LOOP PERTURBATION THEORY 

The goal of our investigations is to see, how the different levels of resummations improve the result. Thus first we 
start with the one loop perturbation theory. Here we need the renormalized Lagrangian; in Feynman gauge it reads 
(using the fact that the photon self-energy is zero): 

C = -^{df,A^f + ^{ido - m)^' - e§^o* + SZmdo'b - SZ^m^^ - Se^Ao^. (17) 
For the fermion self-energy the one loop diagram is the bubble with the contribution: 

- i^Uoo^iP, m) = i-iefj^, iGooik) iGip - k) = -e'u^ j k-\ ^e p, - k,^- m ^ ie ^'^^ 
Moreover we have wave function and mass renormalization counterterms 

5^cf (p) = -SZpo + dZm-m. (19) 
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The complete one loop self-energy is Ei/oop + 5^ct- In the calculation we have to take care of the non-standard form 
of the free fermion propagator. The details of the computation can be found in the Appendix, as a result we obtain 



T,iioopip, m) = —{po - m) 

TT 



-In 



where a = e^/(47r) and 



2?, -l + l + i(ln7r-7£). 



(20) 



(21) 



For renormalization we have to subtract the divergences with help of the counterterms, the finite parts are fixed by 
the renormalization scheme. In the MS scheme we choose the counterterms like 



this results in 



AV) = (Po - m) In . 

TT [l 



The discontinuity of the renormalized self-energy reads 



Disc S(p) = 2a (po — w)6(po — Ti). 

Po 



1 



1 



For the one-loop propagator we obtain 

Q{v) = 



This is consistent with the exact result ([T4l) in the leading order of e} . 
The spectral function DisCp^ iQ reads 



TT [l 



e{v) 



Q{po ~ to) 



2a 



a po-m 
1-1 in 



a 



This spectral function is normalizable, since 



oo 

f dpo TT 

J 2tt a 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



On the other hand the one-loop result is not reliable when | ln(po ~ ^ ic- in the vicinity of the mass 

shell as well as in the large po regime. In order to have a better description of these kinematical regimes, we need 
resummation of certain class of diagrams. 



IV. 2PI RESUMMATION 



As it is discussed in the Introduction, the next level of our approximations is the 2PI resummation. The idea is to 
use the exact propagators in the perturbation theory, this propagator is determined self-consistently using skeleton 
diagrams as resummation patterns. The one-loop bubble diagram in the present case generates the resummation of 
all the "rainbow" diagrams. To obtain an expression for the 2PI resummation we use the technique of [6*1 : we use the 
Hoop formula (jl8l) . interpret the appearing propagators as full propagators, and finally perform renormalization with 
the same form of divergent parts of the counterterms as in the 1-loop case (the actual values will be different). 

The tree level photon propagator is exact, therefore we can write 

^ f d*k g(p~k) 
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Using a spectral representation for the fermion propagator (using that now the Feynman propagator is the retarded 
one and that the fermion spectral function is g{uj < 0) = 0) we find 

oo 

Hp) = -«e J —J k^+iep„-ko-uj + ie- ^ ' 



From this form it is clear that we obtain the weighted one-loop result, ie. 

oo 

— g{uj)Euoopip,uj)- (30) 



In particular, if g{u!) = 2'k5{ijJ — m), then we get back the one-loop result. 

At this point it is worth to examine the UV divergence structure of the 2PI approximation. UV divergences may 
occur in ([30|) for large values of w: using (|20|) we find that the large lo behavior of the one-loop self-energy reads: 

^iioop{p, ^) = -u^ fin --v\+-(- In - + 2?,) po + O(^). (31) 

Since g is integrable for large uj values, therefore the 0{ui^^) is already finite. Therefore the divergence structure of 
the self-energy is A -t- Bp^, just like for the free case, and so the same type of counterterms are needed (although the 
values are different). This is a manifestation of the general case of counterterm renormalizability of 2PI resummations 
f6|. 



A. Analytic study of the 2PI equations 

First we try to analyze ((30)) with analytic methods. We differentiate it with respect to po to find 

dT.uoop a( oj-pQ-ie \ d'^^uoop a 1 9^1; a 

—5 — = - -In 1 , — r~' ^7^ = ^- 

opo TT \ fi J apQ TT po — u! + ie ap^ tt 

Since = pq — m — Y,, we find for Q^^: 



dpi 



= -. (33) 



To solve the equation we first should realize that the (2 = and a ^ cases are very different. If a = then 
^^-ly/ _ Q g^j^^ ^]_jg propagator behaves as ^ = Z /[p^ — m) with some wave function renormalization constant Z and 
mass TO. This agrees with the free case. We also see that the integration constants correspond to the renormalization 
scheme (here the wave function and mass renormalization). 

If a 7^ then we can redefine the variables with an arbitrary Qo scale as 

Po), 1' = -aoa"\ (34) 

V Vl 

then we find 

This equation does not depend on the coupling any more. The coupling constant dependence shows up in the 
integration constants which are the manifestation of the renormalization scheme. We shall also note that the equation 
does not give information about the sign of E and 5*, because for E — ^ —E or — ;> the equation remains the 
same. The chosen signs in turn out later to be the physical choice. 
We introduce 

dy dy d'^ dy , , 

y^iE ^ lE^d^iE^^'d^- 
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This means that we can write for y: 



2yvl/J| = l ^ y=^ = Vhr¥ + yo, (37) 



with an integration constant y^. Therefore 



Vln^" 



= E. (38) 



There could appear an integration constant also here on the right hand side: E — Eq. But recalling that E ^ pq ~ m, 
we see that Eq corresponds to a mass shift: if the mass remains the tree level m then Eq = 0. 

This is the (implicit) solution of the 2PI equations. We see that for real 4* the left hand side is real and positive, 
moreover for ^'(i? = 0) = 1. The E < part corresponds to imaginary values of Vl/. Since the equation itself is real, 
if ^E* is a solution, it is too. This means that the imaginary part is in fact the (half) discontinuity of the solution. 

We see that irrespective of the value of yo, at E — 0, ie. on the mass shell = 1 and so Q = —Qq finite. This 
yields difficulties when we try to apply renormalization conditions on the self-energy. Namely, if we keep the mass 
shell unchanged (this would correspond to the choice of Eq above) , then the renormalization of the self-energy would 
mean S(po = 'm-) — and S'(po — 'm-) — finite. Then, however, near the mass shell the propagator should always 
behave as l/(po ~ "n^): ie. infinite at the mass shell. This means that the physical renormalization process requires 
Qo ^ oo. In this case the propagator behaves near the mass shell as: 



^ 2a, , Po-m 

i + Qoyoy —(m-po) 

because if 4' is close to 1 then the log term can be neglected in p8l) . and we find = 1 + yoE. 

For large values of ^E*, on the other hand, yo can be neglected. Then the integral can be evaluated as 

0Ferfi(\/hr*) = E. (40) 

For large ^ values it behaves as 

^= = E, for large (41) 



Numerical solution 



Now let us turn to the numerical study of the system, based on [6| and [2^ : we determine the discontinuity of the 
self-energy self-consistently. The discontinuity of ((30|) now reads 



(42) 



po 

a f 

DiscE(p) — — duj{po —u})q{uj). 
Po n J 



Knowing the discontinuity of the self-energy, we can use the Kramers-Kronig relation to restore the complete self- 
energy: 



I duj Disc^ i^{uj, k) 

^(P) = / o 7~- — ■ 

Ztt Po — uj + i£ 



(43) 



While is a completely finite expression, in the Kramers-Kronig relation we will find divergences. This corresponds 
to the divergences of the self-energies which must be made finite by applying the appropriate counterterms. Technically 
one can regularize the integral in (j43p and then make it finite with counterterms, or use the (twice) subtracted form 
of the Kramers Kronig relation. To see how it works, we determine the one-loop result from the_tree level spectral 
function and the dimensional regularization of the Kramers-Kronig equations (interpreting u - 



2a/i 



2s 



Jl-2e, 



(2^)1- 



■2e 



uj{po-m-uj) ^ = — (po - m 

ZTT 



2 in h In TT 

e fj. 



(44) 
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The divergence structure is the same, and also the MS scheme result is the same as in ([20)1 (the different finite parts 
are due to the different regularization method). 

Now we can set up an algorithm to solve ([^^ . We choose an arbitrary spectral function as a starting one (practically 
the free spectral function), then follow the following steps: 



step 1:: compute the discontinuity of the self-energy using (|42|) 

step 2:: compute the complete self-energy using the Kramers-Kronig relation (|43 



step 3:: renormalize the self-energy with local counterterms. To fix the counterterms we used on-mass-shell (OM) 
renormalization scheme, ie. the real part of the self-energy at the mass shell is zero and its derivative is also 
zero 



TJ \ a rfReS(po) 

Ke Zj(po = m) = 0, 



dpo 



= 0. (45) 

po=m 



We note here that releasing the first condition yields a mass shift, releasing the second condition yields a 
finite wave function renormalization. But in all renormalization schemes it will remain true that near the 
(renormalized) mass shell the propagator behaves as G{pa ~ m) — C/iPo ~ "^)- 

step 4:: construct the new spectral function from the discontinuity of the propagator knowing the real and imag- 
inary part of the self-energy as 

2ImS(p) 

" {po-m - RcI](p))2 + (ImS(p))2- ^^^^ 



step 5:: continue with step 1 until the process converges. 

Integrations in the above algorithm are performed numerically. This strategy was applied successfully for the $^ 
model in Q. 

The direct application of this strategy, however, this times fails. Numerically what we can observe is that the 
spectral function becomes more and more shallow, and pointwise it goes to zero lim„ Qn{p) ~ 0. In order to see a 
convergence, we had to use a supplementary step in the iteration after step 4: 

step 4':: use a rescaling of the generated spectral function: 

q{p) ^ Ag{Bp) (47) 

with appropriate A and B which can ensure convergence. 

The appropriate values can be found by inspection, but the actual values are not too important (we used A = 73 and 
5 = 11 in our numerics). In this way finally we succeeded to see convergence in the spectral function. 

The numerical reason of this behavior is that the exact spectral function has a discontinuity at the mass shell, and 
- apart from this single point - it has always negative derivative. Numerically, however, we cannot have a jump, 
since in all rcgularizations equation (|42l) yields g{po « to) ~ {po — to)" where n > 2. Since the exact curve starts 
to bend downwards, the recursion tries to lower the spectral function in order to have smaller derivative near the 
mass shell. Since the spectral function has to be positive, these requirements can be satisfied only with g = 0. With 
the continuous rescaling we can achieve that the numerically badly conditioned part, the vicinity of the mass shell, 
becomes smaller and smaller. 

The numerical results can be seen on Figure [TJ The expected asymptotics can be nicely identified on the calculation 
(cf. Fig. [5]). This also proves implicitly that the strategy to resolve the aforementioned numerical problem with the 
2PI equation was correct. 

If we compare the 1-loop, the 2PI and the exact results we see that there is not too much improvement. The IR 
problem near the mass shell which made the 1-loop calculation unreliable, seems to be cured, but in fact the result is 
not closer to the exact one as the one- loop result. The physics of the deep infrared photons cannot be described by 
the 2PI approximation. 
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FIG. 1. The real part and discontinuity of the 2PI propagator. 
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FIG. 2. The expected asymptotics plotted on the data 



V. SCHWINGER-DYSON EQUATIONS AND WARD-IDENTITIES 



The next level of the approximations is based on the Schwinger-Dyson equations. For the Bloch-Nordsieck model 
in Feynman gauge it can be written as 



(48) 



where is the vertex function. 

For the vertex function there is another exact equation, coming from the current conservation. This results in the 
Ward-identity analogous to the QED case [23 |: 



k^T''{k;p-k,p) - g-\p) - g-\p - k). 



(49) 



In this model, however, the vertex function is proportional to . In principle the Lorentz-index in this model can 
come from or from any of the momenta. But, since the fermion propagator depends on the 4- momentum in the 
form UfjtP^, the fermion-photon vertex does not depend on the momentum components which are orthogonal to u^. 
Therefore the Lorentz-index which comes from in fact comes from the longitudinal part of g'*, ie. proportional to 
M^. So we can write r^(fc;p, q) = u^r{k]p, q). 

This gives us the possibility that from the Ward identities we exactly determine the vertex function. The Ward- 
identity for the current conservation yields then in case when u — (1, 0, 0, 0): 



k^Tf'ik; p-k,p)^ koT{k; p-k,p)^ G^^p) - g-\p - k) 



r{k;p-k,p) 



g-Hp) - o~Hp - k) 

ko 



(50) 
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Therefore we find 

m = --^/ ^5{p - k) {g-\p) - g-\p - k)) . (51) 

This is an exact equation in the Bloch-Nordsieck model. Now we will solve this equation in the renormalized theory, 
and demonstrate that the solution is indeed identical with the Bolch-Nordsieck solution presented in Section [III 
In the second term Q^^{p — k) drops out, resulting in an integral 

because of ko — ^ —ko symmetry. What remains is 

= g-\p)He') J ^ ^Gip - k). (53) 

This form is true in the original model, we shall now find the renormalized form. First we adapt the wave function 
renormalization for the fermionic fields which changes the bare propagator to 1/ {ZpQ — {m-\-5m)) where Z = l+SZ. We 
will assume that the mass shell remains the same, then m + Sm = Zm, and the free propagator will be 1/{Z(pq — m)). 
We will use also the notation eb = e + 5e. The full propagator then reads 

g-\p) = Z{po-7n)~^(p), (54) 

Using ([55)) we find the equation 

G(P) - (55) 

Po — m + te 

where 

CiPo)-'-^^ and JiPo)--^elJ^^^gip-k). (56) 

Cipo) can be interpreted as a running wave function renormalization constant. 
With a spectral representation 



J{po)-J^,{u^)h{^-Po + ^e), where = -^/(|^ ^73^1^ ^ (57) 



In the Appendix we evaluate Ii{pq), and we find 



(58) 



— OO 

We rewrite it into ([55)). then, assuming normalizable spectral function, after some algebraic manipulation we find 

OO 

1 1 ^ 1 f duj uj-po + ie 

h -X>c / — g{uj) In 

aft TT IT J It: fi 

CiPo) = . (59) 

Z/ab 

We may assume that the explicit integral is not UV divergent (it can be checked a posteriori, or, as in the present 
case, knowing the exact solution). Then the above equation can be made finite by requiring 

1 1-,. 1 Z Zr 

— + -I?e = — , — = — . (60) 

where a^. and Zr are finite. This form can be interpreted physically as the appearance of the renormalized coupling 
and the finite wave function renormalization Zj.. We note that the coupling constant renormalization equation agrees 
with the nonperturbative coupling constant renormalization in the 0(N) models j23| . 
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Now we find 



oo 



C(.o)4 1-f f'^,i.)ln^^^^^\. (61) 

\ — CXj 

This function depends on the arbitrary scale /i, but the physics, of course, must be fi independent. This can be 

achieved by appropriately changing the Zr and constant when we change fi. The /i-independence of Cipo) requires 
(using the sum rule for g): 



oo 



dCiPo) 1 dzr ar f duj , . uj - Po - ie \ 1 dar f duj uj - po - is I Ur 

— ' ^ ' g(oj) in — / — g[uj) In 1 =0. (62) 



dln/x z^dln/i I n J 2Tr fi j z^Trdln/i 7 27r /i Zr tt 

This can be satisfied if 

1 dZr ^ 1 _ Q 1 dZr Ur 1 dttr _ ^ ^^^^ 

z^dlnn Zr TT ' z^dln/i TT z^Trdln/i 
The second equation means Zr — a^/ao where ao is a constant; the first equation then reads 

dlnzr Ur dur al 1,1 1 1 M ^ t \ arino) , . 

— = — =^ — j = — ^ ^ H — r = - In ^ ar{^i) = 7 — r (64) 

aln/i TT dluyu tt ar{p-) Q;,.(/Lto) tt ^ "r(Mo) Mo 

TT ^ 

Using the normalizability of g we finally find 

oo 

Cbo) = ^ g{^) In A = ^e^. (65) 

— oo 

The ao and the scale A are renormalization group independent quantities (ie. independent of the scale n), these char- 
acterize the renormalization scheme. The appearance of a scale A is the manifestation of dimensional transmutation. 
Now, instead of that scale A it is worth to use M for which KeC{M) = 0. Clearly M ^ A ii A ^ m. Then with 
differentiating C with respect to po we find 



dCiPo) "0 f duj g{uj) 



M 



dpo Ti" J 2n Po — u) + ie 

— oo 

This gives finally 



-Gipo) ^ C(Po) = ^ J dujg{u:). (66) 



M 



{po - m)g{p) = — I dwQ{uj). (67) 



TT 

PO 



By differentiation with respect to po we find 

{po-m)g' + g = -^g => ^(p) =5o(Po -m)-i-^, (68) 

TT 

where go is an arbitrary constant. This is indeed the solution of Bloch and Nordsieck ([Ti)) . now in terms of the 
renormalized quantities. 

But we also see that the condition g{pQ — M) — can be satisfied only with go — 0. This is in close relation with 
the fact that at the mass shell po ~ rn, the propagator (and its discontinuity) is not integrable. 

The lesson of this analysis is that the deep IR physics is well describable by the Schwinger-Dyson equation, truncated 
in a way which respects the Ward-identities. As we have seen, this strategy is renormalizable and exact in case of the 
Bloch-Nordsieck model. 

A big advantage of this approach is that, besides being exact in the IR, it can be easily generalized to other theories. 
So we expect that in QED the Schwinger-Dyson equations truncated in the way we have done it in the Bloch-Nordsieck 
theory will represent the exact result well in the problematic deep IR regime. 
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VI. CONCLUSIONS 



In this paper we examined the exactly solvable Bloch-Nordsieck model from the poin t of view of different pertur- 
bative methods. We first reviewed the known method to obtain the exact solution [l3|. Then the different levels of 
approximations, like the one-loop level perturbation theory, the 2PI resummation and the truncated Schwinger-Dyson 
equations were studied. The 1-loop result exhibits an IR sensitivity when we approach the mass shell which renders 
the theory ill-defined. The self-energy (2PI) resummation reorganizes the perturbative series in a way that this IR 
problem disappears. This does not mean, however, that the result itself would be closer to the exact one, only the 
explicit IR sensitivity cannot be seen. On the other hand, the Schwinger-Dyson equations, truncated in a way that 
the Ward-identities are satisfied yield the exact result in the Bloch-Nordsieck model. This is a new way of obtaining 
the exact solution in the Bloch-Nordsieck model. And, while the original solution method is very hard to generalize 
to other theories, the generalization of the specially truncated Schwinger-Dyson equations is straightforward. 



ACKNOWLEDGMENTS 



The authors thank useful discussions with T.S. Biro, F. Csikor and A. Patkos. This work is supported by the 
Hungarian Research Fund (OTKA) under contract No. K68108. 



Appendix A: Details of the one loop calculation 



The one-loop contribution to the self-energy reads, with a generic u vector in Feynman gauge: 

d^k 1 1 



• 2 2 

-le u 



(Al) 



This is Lorentz-invariant, if we a Lorentz transformation both on u and p. So we may choose a special frame where 
Au = (uo, 0,0,0). If u is a proper 4-velocity, then uq — 1; if it is u = (l,v), then uq = Vl ~ v^, but still constant, 
since u is a parameter of the theory. We find then 



E = e^uo/o( Po-ie), /o(a) = i 



Mo 



d^k 



(27r)4 k^+i£ a + ko' 



(A2) 



Thus it is enough to consider Iq only. There we transform to positive frequency integrals 



/o(a) = i 



d'^k 



(27r)4 k'^ + ie a + ko tt 



— I dko 



(27r)3 k^-k^ + ie - k^ tt 



- / dko 



d-^k 



(A3) 

where in the last step we performed Wick rotation (the choice of the imaginary part of a is crucial for the direction 
of the rotation on the complex plane). 

Now we can write up the integral in fcp f^n.d k space, in the latter using 3 — 2e dimensions: 



/o = a ^2 



dko f d^-^^k 1 



1 



TT J (27r)3-2e k2 _^ k2 ^2 + ^2 ■ 



(A4) 



We use the relation 
d'^-^^k 



2e 



(27r)'^-2'^ 
to proceed as 

oo 

Io= - / dko 



1 



2(47r^ 



2\e 



(4^)'^/2r(d/2 - e) 



dk k" 



-l+2e 



(4V) 



2\e 



(47r)'^/2r(d/2 . 



dz 



-l-e 



/(^) 



(4^A^2) 



2\e 



+ kl (47r)3/2r(| - e) 



ar(-i + e)r(l-e) M^x^ 



,,,f-i-e(,2_,^).r^«(4V)^r(-i + s) 



167r2-v/7r 



87r2 0F 



i - 2 In - + 2 + In TT - 7_B 



dko^, 



(A5) 



(A6) 
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We write it as 



47r2 



- In ■ 



where 



Therefore 



^ 1 , In TT — 7e 

= h IH ^. 

^ 2£ 2 



o2 r 



2 In h 2 + In TT — 7b 



We also need to compute 
r , . ■ . d^k 1 1 1 



dko 



(27r)4 k^ + ie ko a + ko n J "J {2tt)^ fcg - k2 + is a? - kl a 



d^k 1 



1 1 



(A7) 



(A8) 



(A9) 



47r2 



- In ■ 



(AlO) 
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